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ABSTRACT 

We consider the effects of graviton multiplet fields on transitions between string gas phases. 
Focusing on the dilaton field, we show that it may obstruct transitions between different 
thermodynamic phases of the string gas, because the sign of its dimensionally reduced, 
T-duality invariant, part is conserved when the energy density of the universe is positive. 
Thus, many interesting solutions for which this sign is positive end up in a future curvature 
singularity. Because of this, some of the thermodynamic phases of the usual gravitating string 
gases behave like superselection sectors. For example, a past-regular Hagedorn phase and an 
expanding FRW phase dominated by string momentum modes cannot be smoothly connected 
in the framework of string cosmology with positive sources. The singularity separates them 
like a geometric precipice in the moduli space, preventing the dynamics of the theory from 
bridging across. Sources which simultaneously violate the positivity of energy and NEC could 
modify these conclusions. We provide a quantitative measure of positivity of energy and 
NEC violations that would be necessary for such transitions. These effects must dominate 
the universe at the moment of transition, altering the standard gas pictures. At present, it 
is not known how to construct such sources from first principles in string theory. 
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1 Introduction 



String theory differs dramatically from quantum field theory at short distances and high 
energies. At those scales the extended nature of strings comes into play, and may presumably 
regulate ultraviolet divergences inherent in field theory. Yet, since the string scale is high, 
i s ~ x ^, TeV (with the bound being saturated only in most optimistic cases), it is well neigh 
impossible to test string theory in low energy experiments. On the other hand, successes of 
modern cosmology have led to the picture of an expanding universe, which was much denser 
and hotter in its early age. At such high energies, full blown string dynamics may have taken 
place. Its effects may be important for many reasons. In the limits of string theory which are 
under full computational control, calculability requires the presence of additional dimensions 
of space, which somehow must have ended up hidden from low energy probes in the course 
of cosmic evolution. String theory may also have played a role in controlling and selecting 
cosmological initial conditions and subsequent evolution which produced a large and smooth 
universe, by inflation or otherwise. In the course of this evolution, some signatures of string 
theory may have survived the long era of cosmic dilution, and they could yield new insights 
in the microphysics of our world. Ultimately, one may hope that string cosmology can shed 
new light on cosmological singularities, where field theory formulations of matter and gravity 
fail completely, having no direct means to describe the universe at the densest state. 

For these reasons, significant efforts have been invested in the formulation of string cos- 
mology. Interesting models emerged, which treat string matter in thermal equilibrium [1,2,5], 
evolving under the influence of gravity but with specific stringy phenomena accounted for 
in the description of the string matter. This approximation clearly cannot fully reveal the 
details of the microscopic dynamics. Yet, one may at least explore the effects of long range 
string phenomena on some of the outstanding cosmological problems [6-18]. An intriguing 
idea was that at very high energies, because strings are extended objects, the string matter 
can undergo a Hagedorn transition, where the energy density saturates at a finite value con- 
trolled by the Hagedorn temperature, taming the violent divergences in the field theoretic 
models of matter [5]. If this occurs, then using the network of string dualities one might be 
able to describe the universe beyond the string scale by some dual picture, which admits a 
different, but controllable perturbative theory. 

On the other hand, following the growth of the stringy Hubble volume, one may imagine 
that as the universe cools the Hagedorn phase should connect onto the standard expand- 
ing cosmology describing the evolution of dilute strings. The post-Hagedorn stage should 
initially be dominated by relativistic strings, behaving as a radiation dominated FRW uni- 
verse [5]. Similarly, one can also imagine how the Hagdorn phase could be an end product of 
a collapsing universe controlled by string winding modes which are getting dense and trans- 
muting into the Hagedorn phase. If these transitions can be realized consistently, one could 
explore their implications^ for large universes such as our own. The question is whether such 
transitions are possible. 

^ome novel ideas for generating scale-invariant perturbations of the geometry were proposed [19]- [21], 
based on the proposed transitions between different stages of the string gas. However to date they have not 
succeeded [20]. Perhaps it is interesting to explore them further, at least to verify how robust the connection 
between scale-invariant perturbations and inflation is. 
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Our aim here is to study in more detail the effects of the graviton multiplet fields on the 
transitions between the string gas phases. We will focus on the dilaton field, and show that in 
fact it stands in the way of smoothly connecting different thermodynamic phases of the string 
gas. This is because there is a discrete charge characterizing the dilaton sector: the sign of 
the dimensionally reduced effective dilaton, sgn(0 s ) is a conserved quantity, when the energy 
density of the universe is non-negative [22] . The solutions fall into branches characterized by 
this sign, and have a general feature that the location of the cosmological singularity strictly 
correlates with the sign. On the (+) branch, the singularity is in the future, while on the 
(— ) branch it is in the past. The singularity may be either a pole super-inflation, or a Big 
Crunch, which can be reinterpreted as a pole super-inflation in the T-dual geometry. On 
the other hand, if the effective energy density is allowed to dip below zero temporarily, the 
transitions may be facilitated [23] , and for conventional sources which satisfy the null energy 
condition (NEC) those that can complete are the (— ) — > (+) transitions [22]. Such solutions 
describe the cosmologies which develop future singularities induced by the transitions. The 
transitions in the opposite direction, from (+) to (— ) branch, are impossible unless NEC 
is violated [22,24]. This means, that the universes which start evolving towards a future 
singularity cannot evade this fate as long as NEC holds. 

This implies that some of the proposed connections between the Hagedorn phase and the 
expanding FRW phase dominated by string momentum mode radiation [5] are impossible 
in the conventional framework of string cosmology, since those solutions lie on different 
branches. Instead, the decaying Hagedorn stage ends up as a collapsing FRW cosmology, 
which crunches into the future spacelike singularity. Similarly, we will see that a collapsing 
universe dominated by string winding modes will not condense into the long-lived Hagedorn 
stage evolving towards the attractor, again because its dilaton charge is different. The 
thermodynamic phases of the usual gravitating string gases behave like superselection sectors, 
characterized by whether they evolve towards or away from a singular region, which separates 
them like a geometric precipice in the moduli space, where the dynamics of the theory cannot 
bring them together. 

Sources which simultaneously violate the positivity of energy and NEC could modify these 
conclusions [22,24]. There are ways to mimic violations of positivity of energy in string the- 
ory, such as by inducing negative cosmological constant, allowing for Casimir energy effects 
or dimensionally reducing on compact spaces of positive curvature [24]. To evaluate the pos- 
sibilities and impossibilities of any NEC violations, we retain the perfect fluid description 
and provide a quantitative measure of positivity of energy and NEC violations that could 
facilitate the transitions between different thermodynamic phases of string gas. These ef- 
fects must be dramatic, as they need to dominate the universe at the moment of transition, 
altering the standard gas pictures [25]. Such dramatic NEC violations are not derived from 
first principles at present. Further, generically there are dangers that covariant theories 
which violate positivity of energy and NEC may give rise to ghosts. Perhaps it is possible 
to tame the ghost instabilities if Lorentz symmetry is broken [26] , and induce cosmological 
bounces away from the singularity as in [27, 28] . On the other hand, string theory in its 
standard formulation is known not to have ghosts [29-31]. One must properly account for 
this when attempting to build cosmological scenarios which strive to employ string matter 
phase transitions as a method for generating cosmological perturbations. Before inventing 
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mechanisms that rely on crossing between various phases of string cosmology, one must show 
that the bridges that need to be crossed aren't too far, but are firmly grounded in theoretical 
foundation. 

2 A Review of Gravitating String Ensembles 

Cosmological solutions notoriously involve spacelike singularities, where the curvature blows 
up, and the effective field theory description of geometry breaks down. In string cosmology, 
one hopes that the singularities of the effective description could be ameliorated by string 
effects. The complexities of string theory in time dependent backgrounds are a formidable 
hurdle to this program, having so far obstructed singularity resolutions, even casting doubt 
on their existence. To gain some insight in what may be going on, one may instead seek 
to improve the effective theory description by resorting to some fundamental features of 
string theory, such as dualities, which may even survive supersymmetry breaking. Indeed, 
in the string gas approach one should follow the effective geometry to near the singularity, 
and then invoke a hypothesis that as the original geometry begins to explore curvatures 
beyond the cutoff, one may replace it by a dual geometry, where the curvatures again turn 
up small and a perturbative description is restored [2,5]. This very appealing idea thus offers 
a way to address singularity resolutions entirely within the effective theory language. In this 
framework one can ask questions about the consistency of duality-related solutions with the 
desired asymptotic limits to see if they are at least compatible with the standard lore of 
string theory. 

This is the philosophy which we will pursue in what follows. Assuming that we begin 
in some large universe where effective theory works well, one can coarse-grain the sources 
of gravity that inhabit the geometry, using the rules of string thermodynamics to determine 
the leading-order equations of state and gravitational properties of the string ensembles. 
This yields the evolutionary laws for homogeneous cosmologies, which we can follow to 
near the strong coupling/large curvature regime, where the string apparent horizon becomes 
comparable to the string scale £ s ~ M s _1 . At this stage, the interstitial distance between the 
microscopic strings becomes small and the coarse-grained description of cosmology must fail. 
However, if we follow the dictum of dualities and take it that as distances become smaller, 
there exists a dual description where the relevant degrees of freedom are getting more dilute, 
we can assume that it also admits a coarse-grained description and simply ask the question 
about how these two descriptions connect together. In this way, all we need to do is to 
correctly and systematically follow the (semi-)classical solutions of the theory based on the 
truncation of string effective actions to graviton multiplet two-derivative sector to near the 
string scale [5]. In the string frame the action is given by 

where C m encodes the string ensemble contributions, the effects of stringy corrections, and 
possibly any other sources of energy density, including the dilaton potential. The dimensional 
scale k n = Mg~ N is the string scale in D = N + 1 dimensions, and S is the string frame 
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dilaton, related to the string coupling by g 2 = e^ s > . The subscript s indicates a quantity 
computed in string conformal frame, and not only a string gas source. For the exploration 
of gravitational effects of coarse-grained string sources, it is enough to explore the theory (JTJ) 
only on homogeneous backgrounds. However, to elucidate the form of the string gas source 
terms, following [5] one can use the anisotropic metrics 



N 

v(i) 



ds 2 = -n(t) 2 dt 2 + e 3 {t) , (2) 



i=l 



where ai l \t) = e^'^ are the string frame scale factors of the toroidal directions parameter- 
ized by x l s, and n(t) is the lapse function, which can be gauge fixed to unity. While such 
configurations are generically infested with curvature singularities [5] , where some dimensions 
shrink to zero while others explode, away from the singularities one can use thermodynamic 
arguments to deduce the form of the string sources, and specifically energy density and pres- 
sures, at large scales in the gravitational equations. So, specializing for the moment to only 
ensembles of gravitating string gas, dimensionally reducing the theory (OQ) on the commuting 
spatial Killing vectors d{ of (E]), and using the shifted string frame dilaton <p s = 4> s —^2iL\ As , 
one finds the effective ID action [5], 

^-/*»fe(E^-0)+^".»«}. (3) 



where F is the comoving free energy of the string gas, which depends on the temperature 
f3^ 1 and the scales of the comoving 'container' A s . Varying this action with respect to 
intensive parameters and using the standard relations between energy, pressure and free 
energy, E s = —(F + /3dpF) and P { = d^F = -(d x (i>E)\p =const , in the gauge n = 1 one finds 
the string frame Hamiltonian constraint and equations of motion 

N 2 

AS -^?) = \e*-P®, (5) 

= (6) 



i=i 



2 



respectively, where we have chosen the units M s = 1. The sources obey the conservation 
equation, 



N 



E s + ^A«P, = 0, (7) 



i=i 



which follows from (jl])-©. We will focus on the sources of the form E s = YliLi E s o e/3A ° ' > i n 
the limit when the distribution is isotropic, As = A s . The pressures in different directions 
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Pi are replaced by the mean pressure P s = J2i=i ^» = ]v Si=i = ~n^ s ^ sucn ^ na ^ 

the equation of state parameters 7W = pjf 1 /E s reduce in the isotropic case to 7 = —f3/N. 

Once the equations of motion (HJ)- (J6]) are derived, it is convenient to replace the comoving 
energy E s and comoving pressure P s by the energy density p s and pressure p s . They are 
defined by p s = E s e~ NX and p s = P s e~ NX , such that the relevant sources in the equations of 
motion (JIJ-Q can be substituted according to 

e*'E a = e+'p a% e^P 8 = e*'p B . (8) 

This allows one to rewrite the equations (SJ)— (J6J) in terms of the original N + 1-dimensional 
dilaton variable (fi s as 

(^2^)A s 2 + ^-iVAA = \e+-p. t (9) 

A, o.X I A A; = \e+'p a , (10) 
p s + N\ s (p s + p s ) = 0. (11) 

The last equation is the stress-energy conservation, which follows from ((Tj) and jSJ). The 
dilaton equation of motion follows from stress-energy conservation, as one can verify by tak- 
ing a derivative of (0) and using Eqs. (I9l- ([TT1) to solve for (fi s . This form of the equations 
of motion is convenient to go beyond the gravitating string ensemble approximation and 
introduce more general sources of stress-energy, which may depend on the dilaton field, in- 
cluding any dilaton potential on the landscape. In that case, of course, the simple equation 
of state p s = jp s will be replaced by a more complicated expression, derived from the La- 
grangian. This will be useful to extend our results to more general sources, particularly when 
formulating the conditions for the connection between different gravitating string ensembles. 

2.1 Isotropic Single Fluid Cosmologies 

The equations of motion (JI])-© or (I9l- (TTT1) can be easily solved for the single fluid com- 
ponent, with a constant equation of state 7. These solutions serve as the attractors of 
the homogeneous and isotropic multi-component fluid cosmologies. Indeed, starting with 
a mult i- component source, an expanding universe eventually ends up being dominated by 
the component with the least positive equation of state parameter 7. Conversely, the fluid 
with the most positive 7 will asymptotically dominate in the contracting universe limits. 
When the fluids are composed of unstable states which may decay, the equation of state 
may eventually change, transmogrifying the dominant source of stress-energy and the uni- 
verse it inhabits into a different fluid model. One may therefore try to patch together general 
cosmologies from such attractor geometries. The key for this quilt work is to determine the 
rules for connecting solutions. At the leading order, one may assume that the decay time is 
much faster than the Hubble time at the moment of decay, and approximate the transition 
by an instantaneous one. In this limit one can again ignore the microscopies and retain 
the long range, coarse-grained description, using Israel junction conditions to relate the field 
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variables at the spacelike boundary between different epochs. We will pursue this approach 
in what follows. 

The starting point is to taxonomize the various interesting single component solutions. 
To count them up, one can first rewrite the system of equations fl9l)- f|TTl) in the canonical 
form, by solving the Hamiltonian constraint (Q. This yields 

S = NH S ± y/NH* + e*°p a , (12) 
H s = ^e+- P .±H.y/NH? + efi' Pa , (13) 
Ps = -NH s (p s + Ps ) , (14) 

where we have introduced the notation H s = X s . We immediately note that the solutions 
can be classified by the sign in front of the root, which by the definition of tp s = <p s — NX S 
is precisely the sign of the reduced dilaton flow, i.e. sgn(y? s ), as we can verify by comparing 
dlj) and ([12]) . and by the arrow of time. The sign of Lp s can change only if the discriminant 
NX 2 + e^ s p s vanishes. Clearly, this is impossible if p s > 0. We will discuss this in more 
detail later on. The sign of H s can change if the universe 'bounces', which may occur in the 
string frame due to the fact that the dilaton kinetic term is not canonical. At any rate, this 
discussion shows that for any value of the source terms there are four solutions, described 
by the combinations of the two signs. The latter sign selection is the same one encountered 
in the pre-Big Bang cosmology, and so we adopt the same terminology here, as it obviously 
generalizes. From now on, we will refer to the (±)-branch by referring to the sign in front 
of the square root. 

Equations (fT2"]) - (fI~3"]) show that as long as p s > 0, in the expanding universe on the (+) 
branch, the string coupling g s = e^ 3 '' 2 increases. We also see that on the (— ) branch, the 
string coupling increases during the contracting phase. On the other hand, the evolution 
of the string coupling in a (+) branch contracting universe and a (— ) branch expanding 
universe is more sensitive to the sources, and depends on the time dependence of e^ a p s . In 
particular, it may occur that the p s > fluid can flip the sign of S in these cases, as can be 
seen by rewriting Eq. (fT21 as <p s = ±\/N \J <fi 2 — e^ a p s + (fi s . This also shows that if p s < 0, 
the coupling is fixed to grow on the (+) branch contracting solutions, and decrease on the 
(— ) branch expanding solutions. 

We should note that we can generate new solutions by the symmetry transformations of 
the effective action, which are the time-reversal dt — > —dt and T-duality [4]. In fact, any 
two solutions which reside on the same branch, as defined by the sign of (p, but differ by the 
sign of H s , are T-dual images of each other. The T-duality operates on the background as 
the map A s — > — X s , (p s — > ip s . This also implies the change of sign of the equation of state 
parameter 7 = p/p, which is easy to see from the fact that the T-dual of the decaying energy 
density increases, and vice-versa. An important thing to stress is that T-duality does not 
change the branch since the sign of dt, and therefore (p s , remains unchanged. On the other 
hand, the time-reversal operation changes the branch on which the solution resides, since 
under it the sign of (fi s changes. Thus starting from any one solution, we can get four different 
ones by the sequence of maps T-duality — > time-reversal — > T-duality. In the case 7 = 0, 
or when p s = p s = as in pre-Big Bang cosmology [6], this exhausts all the independent 
families of solutions, but not in general cases, as we will see in what follows. 
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Cosmic expansion and/or acceleration in the string frame do not imply the same in the 
Einstein conformal frame, where the graviton kinetic terms are normalized canonically. On 
the contrary, when the dilaton evolves fast, the picture in the Einstein frame may be com- 
pletely reversed, so that an expanding string frame universe looks like a collapsing universe 
crashing into a singularity. This is relevant to understand the obstructions to pasting to- 
gether different cosmic epochs. To go to the Einstein frame, one performs the conformal 
map to the action where the graviton and dilaton terms are canonically normalized. The 
field redefinition which relates the Einstein frame and the string frame variables is 




which on FRW backgrounds yields Further, one must be careful with the 

comoving time relationship as well, which is given by gauge-fixing the two lapse functions 
to unity. This yields dr e = e~^ 3 '^ N ~ l >dt. Using <fi s = <p s + NX S , the logs of the scale factors 
and their time derivatives are related by 

iv^r ' d^=~ { ^^ ] ~- (16) 

From the second of (|T6|) we see that the sign of the Einstein frame Hubble parameter H e = ^ 
can be reversed relative to the sign of H s , depending on the evolution of <p s . This can be 
seen most easily on a case to case basis. 

Let us now write down the exact single fluid isotropic solution, defined by a constant 
equation of state parameter 7 = P s /E s = p s /p s and the scaling law E = Eoexp(-'yNX). 
For this purpose it is useful to go back to the (isotropic limit of ) the equations dlJ-Q and 
introduce a new time coordinate x defined by E = The equations of motion become 

ip' s 2 -NX' s 2 = E^e^ NX % (17) 
A'; - 7 iVA; 2 - V ' S X' S = ^E-^^, (18) 

tf-lNtfX-NK* = \E^e^ NX % (19) 

respectively, where the prime denotes a derivative with respect to x. The solutions of these 
equations are [6,32] 

X s = \ M + lhi[x(x-x,)] + —=hi(l-—), (20) 



Vs = <Pso - - In [x(x - x # )) - 7 — In ( 1 - — ) , (21) 

a a V x 



h = ^In^g-g.)]- ^ 1} VN In (l - ^) , (22) 

a ' a \ x J 

where a = 1 — N^ 2 . The integration constant A s o is pure gauge, and can be changed 
arbitrarily by a rescaling of the spatial coordinates x k — > (x h . The constant ip s o is linked 
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to the comoving energy E according to e Vs0 = ^ E e~ lNXs0 . Note that in principle, we can 
always write the solution in the form where the integration constant is positive. Indeed, 
to do it, it is sufficient to shift x — > x + x*. In this case, the only change in the solution 
( 120]) - (122]) is the sign flip of the terms oc In (l — It is quite clear that the times x = or 
x = x* are special moments in the evolution of the epoch described by (120]) - (122]) . as there 
the scale factor a = e Xs can vanish or diverge. To see precisely what goes on, however, it is 
more convenient to study separately the cosmologies of different relevant string phases. 

2.2 Hagedorn Phase 

In string theory, the density of states diverges at some high temperature, where the theory 
undergoes Hagedorn transition. Beyond this scale, the injection of additional energy into 
the system does not affect the temperature, which remains constant, but simply increases 
the entropy [33]. Essentially, what happens is that the strings become extremely wiggly as 
they keep absorbing the energy [34]. This saturation of the temperature implies that the 
dependence on the other intensive parameters vanes, and so the pressures vanish. Hence 
the Hagedorn phase is described by the constant comoving energy E = Eq and vanishing 
pressure P = 0, giving the equation of state parameter 7 = [5]. This implies that the 
parameter a as defined in the solution (l20]) - (f22|) is a = 1. Rewriting this solution in terms 
of the comoving string frame time coordinate, by using t = x/E , and setting x*/E = t*, 
we find [5] 

X s = A s0 + -^ln(l-|j , 

<p 3 = cp s0 -ln[t(t-Q] . (23) 

The unshifted dilaton is 




so that the string coupling is g s = ^ (l — j-) 2 ■ 

First we note immediately that the special solution with = describes a static string 
geometry, with A2 = A s o = const and (p s = <p s o — 2lnt, where the runaway dilaton absorbs the 
effects of the Hagedorn source allowing the universe to loiter forever. This special solution 
is not a universal attractor, since it is unstable: while there are configurations which flow 
towards it, representing either universes expanding from zero size or contracting from infinite 
size (as exemplified by the Class III and IV solutions described below), their time- reversals, 
which are perfectly allowed by the time-reversal symmetry of the field equations, describe the 
evolution away from this special geometry, either towards a crunch or an infinite universe. 

Now, the solutions for £* 7^ have classical curvature singularities at t = or t = t*, 
where the scale factor either diverges or vanishes. The interval between and £* is classically 
forbidden, since inside it the logarithms are imaginary. Note that the ordering and does 
not introduce extra solutions because the order can be changed by the shift t — > t + , 
which sets > while flipping the sign of the log term in the equation for the scale factor. 
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This means, that the solutions with t* < are simply T-duals of the solutions with t* > 
for a fixed sign of t. On the other hand, the time-reversed geometries are obtained by 
simultaneously flipping the signs of both t and t* in the solution. The solutions and 
(j24p are perturbatively meaningful descriptions of the following classical string geometries: 

• Class I: The universe expands from finite size, starting with g s <C 1, on the interval 
t < 0, with t* > and (fi s > (evolving away from the t* = repeller) until t = 0, 
where the scale factor and the string coupling diverge, and yield a curvature divergence 
as well. Effective field theory description will cease in the course of evolution towards 
this limit. 

• Class II: T-dual of Class I, which contracts on the interval t < 0, with > and 
(p s > (evolving away from the = repeller) from finite size to a Big Crunch at 
t = 0, while the dilaton flows to weak coupling j s < 1. Geometric description ceases 
in this limit. 

• Class III: The universe slowly expands for t > t# > from zero size to a constant 
size at t — > oo, with ip s < (approaching the i* = attractor) and the dilaton grows, 
S > 0, such that the string coupling becomes strong at late times, and weakly coupled 
perturbative picture must eventually be replaced by a strongly coupled description. 
During the early epoch when a < 1, the geometric description fails. This solution is in 
fact time-reversal of the Class II solution. 

• Class IV: T-dual of Class III, which during times t > £* > contracts from infinite 
size to a constant finite size at t —>■ oo with <fi s < (approaching the £* = attractor) 
and has a finite string coupling, that initially grows (<fi s > 0) but eventually decreases 
(<f) s < 0), evolving to weak coupling regime. Again, the geometric description of this 
solution cannot be trusted beyond a < 1. This solution is time- reversal of the Class I 
solution. 

As explained above, there us no more independent solutions, essentially because X s de- 
pends only on the ratio ^j 2 -. This degeneracy follows from 7 = 0, so that all the gravitating 
Hagedorn phases remain connected by duality and time-reversal. The solutions are summa- 
rized in Table (1). Note that the (+) branch solutions evolve towards the singularity, and 
(— ) branch away from it. The evolution of these solutions is represented in the tp s , H s phase 
plane in Figure (pQ), with the convention that the time flows from left to right. 

The perturbative description is limited by the flow of a solution to short distances and/or 
strong coupling. The geometrical singularities appear particularly dangerous. Of course, 
once we recall the philosophy behind dualities, the singularities may be interpreted as merely 
an indicator of dual branches which the solution must be replaced by, once the scale factor 
dives below the string scale. The question then is if such dual maps relate geometries that 
admit asymptotic description in some supergravity. In particular, one would like to be able 
to have a cosmology whose late time behavior describes a radiation dominated stage, with 
the coupling possibly evolving away from extreme weak coupling, so that the theory may 
yield useful mechanisms for moduli stabilization and asymptote to a phenomenologically 
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Class 


Branch 


Expansion 


Shifted Dilaton 


Time 


Singularity 


I 


(+) 


H s > 


<p s > 


t < 


Future 


II 


(+) 


H s < 


S > 


t < 


Future 


III 


(") 


H s >0 


v? s < 


t > U 


Past 


IV 


(") 


H s < 


0s < 




Past 



Table 1: The four Classes of the Hagedorn solution ([23]). In Class II, if (VN+l)/2 > t/U > 1 
then (f) s > 0, otherwise (j) s < 0. Likewise, in Class III if {\fN + l)/2 < t/t* < 1 then (f) s < 0, 
otherwise 4> s > 0. Class I and II, and Class II and IV, respectively, are related by T-duality, 
and Class I and IV and II and III, respectively, by time-reversal. 




Figure 1: The (ip' s ,H s ) phase diagram of Hagedorn cosmologies. The black dots mark the 
Hagedorn fixed points, (<p' s , H s ) = (±y/~Eo, 0). 
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reasonable regime. To see if such a stage can be preceded by a Hagedorn stage, we must first 
explore the phase space of cosmology of excited dilute strings, to which we turn to next. 

2.3 Dilute String Phase 

Below the Hagedorn temperature, the ensemble of strings behaves differently. If spatial 
directions are compact, then the dilute strings can carry energy in two different channels: 
i) the momentum modes, which describe the propagation of disturbances along the string, 
and ii) the winding modes, which store the energy to stretch the string around the compact 
dimension. In the former case, the momentum mode channels are relativistic at high energies, 
and so their thermal ensemble behaves as a gas of radiation in the string frame, with the 
equation of state parameter obeying 7 = 1/N, yielding E = E$e~ x . The winding modes are 
the T-dual of the momentum modes, and so their equation of state parameter, by the fact 
that under T-duality the pressure changes sign, is 7 = -1/N, such that E = E e x [4,5]. 
The negative pressure reflects the fact that when the dimensions expand, the winding mode 
energy increases since it is proportional to the length of the string. Hence the energy density 
dilutes more slowly: although the number density drops down as is usual, with the one power 
of spatial volume, the energy increases and compensates some of the dilution. Of course, in 
the case when strings gravitate, the whole geometry needs to be consistently subjected to 
the T-duality transformation. 

The homogeneous cosmology dominated by momentum or winding modes is still given 



by the general solution (f20|) - (f22l) but with the parameter a now given by a = Thus 
the solution for the scale factor and the shifted dilaton formally reads 

As = A s0 ± In [x(x - x*)] + In (l - ^) , (25) 

N , , y/N ( x*\ 

(fs = ~ n - 1 ^ X ~ T N - 1 V 1 ~ ~x~) ' ^ 



where the upper sign formally describes a momentum mode dominated universe, and the 
lower sign a winding mode dominated universe. The singular points of the solutions are as 
before and x*, and the interval between them is classically forbidden, by the condition of 
reality of the solution. The unshifted dilaton in a universe dominated by the momentum 
modes, which evolve as string radiation, is then given by 

0- = sO + ViVl n (l-^) , (27) 

whereas in the universe dominated by the winding modes it is 

« = A^T ln W +VF WTT ln ( l - X i)' (28) 

as dictated by T-duality. 

To see that the solutions with the opposite signs of 7 are in fact T-duals of each other, 
we need to shift the 'time' by x — > x + x* and flip the sign of the integration constant x*. 
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With this, it becomes clear that \ s — \ s0 — > — (A s — \ s o) while ip s remains unaffected. In 
fact, this becomes obvious if we write the expressions for the scale factor and the exponent 
of the shifted dilaton, 

a s = e s0 x (Vn-ixvn+1) ( x _ x j (vw-mvw+i) ; 

e V s o/2 

x (VJV-l)(v / iV+l) ( x — x \ v ( v / Jv-i)(v'jv+i) 



W2 



However, unlike before, the solutions depend on different powers of x and X X * SO that 
the sequence of T-duality and time-reversal cannot link all of them directly. Instead, there 
are two different families of linked solutions, which arise because T-duality interchanges 
momentum and winding modes, and flips the sign of 7 7^ 0, breaking the degeneracy found 
in the gravitating Hagedorn gas example. 

To understand the universes which these solutions describe, we need to follow their 
evolution for specific values of x*. As before, there are the special solutions with x* = 0, 
which are 

A s = X s0 ± - ln|xj , ip s = if sQ - ^ N _ l ln|x| , <f) s = <Pso + 2 ^ 1 _ ± 1 ^ ln \ x \ ■ ( 30 ) 

They are clearly T-duals of each other, representing asymptotic attractors for momentum 
mode-dominated and winding mode-dominate cosmologies, respectively. For general 7^ 0, 
there are eight distinct solutions that comprise two families interlinked by T-duality and 
time reversal, which are are conveniently parameterized by the branch sign and the signs of 
7 and x*. They are as follows: 

• Class 1: The momentum mode-dominated universe which expands forever as a power 
law, starting with j, < 1, on the interval x > x* > 0, with (p s < 0, such that 
asymptotically g s — > g s0 = const., evolving towards the x* = momentum mode 
attractor. This is a (— ) branch solution. Near the singularity x = x* the geometric 
description will cease, and the classical geometry there resembles a Big Bang. 

• Class 1 TR : The time-reversal of the Class 1 solution, found by flipping the signs of 
both x and x*. Evolves away from the = momentum mode repeller. 

• Class 2: Another momentum mode-dominated universe, which however resides on 
the interval x < 0, with x* > 0, and which starts out as a collapsing universe but 
undergoes a string frame bounce, and runs off to a super-inflating pole singularity at 
x = 0. This is a (+) branch solution, (p s > 0, starting with some g s0 = const, evolving 
away from the x* = momentum mode repeller towards strong coupling g s ^> 1. 
Near the pole, both the geometric description and perturbation theory in g s break 
down. This solution cannot be related to the Class 1 case neither by T-duality nor 
time-reversal, which is clearly from the existence of the string frame bounce. 



Class 2 TR : The time-reversal of the Class 2 solution. Again, the signs of both x and 
flipped. Evolves towards the x* = momentum mode attractor. 
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• Class 3: T-dual of the Class 1 solution, describing a winding mode-dominated uni- 
verse, which is collapsing from infinite size at x = > to zero as i — > oo, with 
tp s < 0. This is a (— ) branch solution. In the course of evolution, the string coupling 
changes from g s ^> 1 to g s <C 1. Initially, the theory is given by strongly coupled 
winding mode gas, which evolves towards weakly coupled dense gas, described by the 
x* = winding mode attractor, where the geometric picture breaks down. 

• Class 3 TR : The time-reversal of the Class 3 solution. Evolves away from the x* = 
winding mode repeller. 

• Class 4: T-dual of Class 2, a winding mode-dominated universe defined on the domain 
x < (with > 0) which initially slowly expands out of zero size, evolving away from 
the x* = winding mode repeller. It reaches a maximum size due to the resistance 
from the winding modes, and starts to collapse again. In the course of evolution, 
(f s > so this is a (+) branch solution. The string coupling varies from g s <C 1 to 
g s ^> 1, so that near the crunch both the geometric description and g s expansion fail. 

• Class 4 TR : The time-reversal of the Class 4 solution. Evolves towards the x* = 
winding mode attractor. 

Rewriting these solutions in terms of the comoving time t is cumbersome. However in 
the asymptotic limits \x\ — > ±oo it is straightforward to extract the leading order form of 
the solution as a function of t. In fact these limits are exactly the attractors fl30|) . Using 
dx/dt = E e~ Xa , and the fact that f[2"5l) converges to X s — ► ±;^ry hix 2 when |se| — ► oo, we find 
that | a; | ~ |t|( 7V - 1 )/( iV + 1 ) when momentum modes dominate. Similarly, when winding modes 
dominate, \x\ ~ ^(N- 1 )/^- 3 ) for N ^ 3, and \x\ ~ e^ '*' for iV = 3, as is straightforward to 
see. Therefore the momentum mode solutions in the limit x — ► ±oo approach 

2 2 TV 

a s -> a s0 |t|^+ T , -» y?s0 - — - - In \t\ , <p s 4> s0 , (31) 

which is precisely the form of an FRW radiation-dominated cosmology in N + 1 dimensions. 
The fact that the string dilaton <p s approaches a constant in this limit despite the absence of 
a stabilizing potential is a well known consequence of the conformal symmetry of the sources 
which has this attractor behavior [35]. On the other hand, as x — > ±oo the winding mode 
solutions converge to 

a s0 2N , , , AN , . . . . 

fs ~^ fso 77 — o m H> Vs -> s o - -T7 — o In |*| , (32) 



for iV 7^ 3, which changes to 

a s -> a s0 e" £o|t| , -> V 9 ^ - 3-E In |t| , S -> s o - 6E |t| , (33) 

when N = 3. We therefore see that the winding mode geometries describe an inflating or 
deflating asymptotic cosmology in the string frame, depending on the direction of the flow 
of t. 
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Figure 2: The (ip' s , H s ) phase diagram of momentum mode-dominated cosmologies. The 
black dots mark the momentum mode fixed point cosmologies, which in the r variable are 
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Figure 3: The (<p' s , H s ) phase diagram of winding mode-dominated cosmologies. The black 
dots mark the winding mode fixed points, ((p' s ,H s ) = (±^3E /2), ±^y E /Q). 
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These solutions are depicted in Figures ([2]) (momentum mode-dominated cosmologies) 
and ([3]) (winding mode-dominated cosmologies), where we use the time variable r, defined 
by the relation dr = exp( Vs ~'^ NXs ) dt (see below), and assume that the time flows again from 
left to right. The two families of solutions mutually linked by time-reversal and T-duality are 
also summarized in Tables (2) and (3). Again, we see that the (+) branch solutions run to a 
singularity, while the (— ) branch ones evolve away from it. We can now turn to determining 
how and when these solutions may be matched onto the Hagedorn gas cosmologies of the 
previous section. 



Class 


Branch 


Expansion 


Shifted Dilaton 


Time 


Singularity 


1 


(-) 


H s >0 


y? s < 


x > x* > 


Past 




(+) 


H s < 


S > 


x < < 


Future 


3 


(") 


H s <0 


S < 


x > x* > 


Past 


3 TR 


(+) 


H s >0 


<p s > 


x < x* < 


Future 



Table 2: Smooth family, without string frame bounce. Classes 1 and 3 are T-duals of 
each other, as are their time-reversed solutions, describing momentum mode-dominated and 
winding mode-dominated cosmologies, respectively. 



Class 


Branch 


Expansion 


Shifted Dilaton 


Time 


Singularity 


2 


(+) 


H s <0- 


*H S >0 


ip s > 


X 


< 0, x* 


> 


Future 


2 tr 


(") 


H s >0- 


^H s <0 


0s < 


X 


> 0, x* 


< 


Past 


4 


(+) 


H s >0- 


*H S <0 


ip s > 


X 


< 0, x„ 


> 


Future 


4 TR 


(") 


H s <0- 


+ H S >0 


v? s < 


X 


> 0, 


< 


Past 



Table 3: A family with the string frame bounce. Classes 2 and 4 are T-duals of each other, as 
are their time-reversed solutions, again describing momentum mode-dominated and winding 
mode-dominated cosmologies, respectively. 



3 Branch Changes and Energy Conditions 

The transitions between different phases of string are normally precipitated by the change in 
internal energy. As the dilute string gas is heated to the Hagedorn temperature, it undergoes 
the Hagedorn transition, and starts to absorb energy by increasing its entropy while keeping a 
constant temperature. Conversely, if the internal energy is dissipated, by decreasing entropy, 
the Hagedorn gas will reach the transition point and begin to cool. An obvious arena for such 
transitions is an expanding universe, where the expansion dials the internal energy of the 
matter up and down. Thus in string theory, it is expected that in cosmological backgrounds 
the gas of string will scan through various allowed thermodynamic phases as the universe 
evolves. This has led [2, 5] to propose models of early universe which may even provide 
a way around the cosmological singularity, not directly by a bounce but rather by a dual 
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relationship between collapsing and expanding universes. In these works, typically one seeks 
to have cosmologies which look like a patchwork of different thermodynamic phases of the 
string, that evolve into each other dynamically. 

In the presence of string gravity, however, the transitions between different string phases 
become subtler, complicated by the presence of the dilaton. The dilaton may significantly 
affect the evolution of the background geometry inhabited by the string gas. If one treats 
it too cavalierly, one may not notice that the dilaton can obstruct transitions between some 
of the gravitating thermodynamic string ensembles. One may argue that just before such 
transitions the dilaton is stabilized by dynamical effects [14,19]. But if the theory is to retain 
its duality symmetries at the fundamental level, the dilaton stabilization and decoupling 
should only occur at low energies, below which the dualities are obscured by the infra- 
red effects. Hence one expects that any phase transitions at high energies, in a very early 
universe, should be subjected to the selection rules dictated by the full, T-duality symmetric 
dilaton gravity which governs the geometry. Let us now determine these rules in a precise 
manner. We will show that as long as the energy density inhabiting the universe is non- 
negative, branch changing is impossible, and the transitions can only occur between string 
phases on the same geometric branch. Further, if positivity of energy is not maintained, 
branches may change, but unless the sources also violate NEC, which in homogeneous and 
isotropic cosmologies requires p + p < 0, it is impossible to go from a (+) branch solution 
to a (— ) branch one. Since (+) solutions have future singularities, this result is really a 
singularity theorem in a different guise: it implies that the solutions cannot be relieved of 
strong curvatures unless NEC is violated. 

Indeed, when discussing the equations (fT2l) - (fT4l . we already observed the first part of this 
statement, that the solutions on opposite branches, as characterized by sgn(0), cannot link 
up unless the discriminant N\ 2 S + e^ a p s vanishes. Since the first term is positive definite, this 
implies that if the energy density p s is also positive definite, a transition between different 
branches is impossible. Hence to have any chance for changing branches, the energy density 
must dip below zero, which of course by itself still does not guarantee that the branch change 
will actually occur. 

Although very simple, this statement of positivity of energy excluding branch changing 
is already quite strong: it already excludes some of the phase transitions described in the 
literature. For example, [5] describe a universe which expands away from the self-dual 
point on the Hagedorn phase, and eventually below the Hagedorn transition it transits to a 
radiation phase, or possibly to an expanding winding mode phase. Near the self-dual point, 
the universe just hovers, with the background asymptotically evolving out of the attractor 
solution, given by the = limit of (123]) . Prior to the self-dual point, the universe may 
have evolved as a time-reversed cosmology of either of the two post self-dual point options. 
However, the inspection of the Hagedorn phase solutions, in Table (1), immediately shows 
that the Hagedorn cosmology that slowly flows out of the fixed point at the self-dual radius 
is the Class I geometry, which is a (+) branch solution. On the other hand, the radiation 
cosmology which has no future singularities but evolves towards a conventional FRW universe 
with a frozen dilaton is the Class 1 geometry, a (— ) branch solution, as we see from Table 
(2). Thus if the energy density dominating the universe is always positive definite, the decay 
of the Hagedorn gas in a cosmology evolving from the attractor will not yield an asymptotic 
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radiation-dominated FRW universe. Alternatively, the Hagedorn gas stage is really evolving 
towards the strong curvature regime, whereas the radiation FRW universe is not. To change 
from one to another, the evolution to large curvature must be averted, and that is impossible 
with the positive definite energy density. 

If we turn to winding mode cosmologies, we can see by a similar argument as above that 
the Class I Hagedorn gas universe, which flows out of the attractor, can only connect to 
the Class 4 geometry, which is also a (+) branch solution (see Table (3)). This transition is 
possible early on while the winding mode universe is still expanding. However, the expansion 
ceases in this universe and it collapses back to a future singularity. Once it passes below 
the self-dual radius, we may replace it by its T-dual description, however this is the Class 
2 momentum-mode dominated cosmology which experiences a pole singularity in the future 
as opposed to asymptoting to the conventional radiation-dominated FRW universe. The 
bottomline is that the strong curvature singularity in the future is unavoidable. Similar 
obstructions can be inferred for the time-reversed configurations describing the universe 
before the self-dual point. 

Thus as long as the energy density is non-negative, we can only sew together the solutions 
on the same branch [22]. For instance, requiring that the late universe is a radiation- 
dominated FRW universe of Class 1 confines one to the (— ) branch solutions. If the evolution 
is then followed backward in time, it is plausible that this universe can be preceded by the 
Hagedorn cosmology given by the Class III solution, which is on the (— ) branch too. However, 
in the far past, this solution does not asymptote to the Hagedorn attractor configuration, 
but instead has a Big Big singularity in the string frame. It is conceivable that the small 
universe region near this singularity can be resolved as the T-dual geometry of Class IV, 
which in turn could have evolved from the Class 3 collapsing winding mode-dominated 
world. In this chain of events, one might be able to 'relegate' the problem of cosmological 
singularity, reinterpreting it as the problem of finding the correct dual description of a small 
universe. On the other hand, in contrast to the discussion of the cosmic phase quilt that 
was discussed in [5] , the universe of this example does not spend an infinite time lingering 
around the Hagedorn attractor. In fact, the Hagedorn attractors are completely excised, 
now being in the future of the transition to FRW radiation, or in the past of the winding 
mode condensation. One could now imagine that the local perturbations in the collapsing 
T-dual geometry, if such is consistently constructed, can be computed, and that they could 
map onto long distance perturbations in the small universe near the singularity. This would 
resolve their short distance dynamics. However this would be completely different from the 
ideas for generation of perturbations in string cosmology which have been pursued so far. 
Nevertheless any applications of such ideas to phenomenology must await the development of 
the framework that might enable us to compute how the local perturbations in one universe 
T-dualize to large scale perturbations in another. The present techniques only allow a 
precise construction of T-dual geometries when there are isometries, which are of course all 
generically broken by the perturbations. 

To understand the dynamics of the solutions with non-negative energy density, it is 
convenient to go to the phase space description. The phase space evolution of solutions can 
be illustrated in simple cases of Hagedorn, momentum mode or winding mode universes, with 
7 = 0, ±l/iV, respectively, in Figures (HI) -(IB]), where we now also plot the flow fields which 
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Figure 4: The ((p' s , H s ) phase diagram of Hagedorn cosmologies, with phase space flow and 
the limiting envelopes (dotted lines) describing the case E — 0. 




Figure 5: The ((p f s , H s ) phase diagram of momentum mode cosmologies, with phase space 
flow and the limiting envelopes (dotted lines) describing the case Eq = 0. 
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Figure 6: The ((p' s , H s ) phase diagram of winding mode cosmologies, with phase space flow 
and the limiting envelopes (dotted lines) describing the case E — 0. 

indicate the direction of evolution for nearby trajectories. In the plots, we start with the 
Eqs. ( jT2l) -( fr4"|) . rescale the time variable to dr = e <y '' f ' s ~'~ lNXs ^ 2 dt and introduce the variables 
I = A' and / = ip' to recast these equations as a first-order system, 

f ~ \f{lNl + f) = ~Eo , V - h(jNl + f) = lEo , f - Nl 2 = E . (34) 
This system has fixed points 

(/•') = (^^^N-^J^k-y) • < 35 > 

which are precisely the relevant attractor/repeller solutions that we discussed previously. 
The flow fields are sensitive to the location of the fixed points, as controlled by the initial 
comoving energy Eq. However, to understand the dynamics we must keep in mind that the 
flows are constrained by the last of Eqs. ( 1341) to follow a specific trajectory in phase space 
as proscribed by the Hamiltonian constrain from the gravitational sector, which in simpler 
parlance is just the Friedman equation. An arbitrary trajectory in phase space for a fixed 
energy density and pressure, that may be a solution of the first two of (IMl) is not a solution 
of the full gravity equations unless the constraint is satisfied explicitly. So, the flow fields 
in the Figures (HI)- ([6]) are to be understood merely as an indicator for how the system will 
evolve if it is perturbed by a source that is asymptotically subleading to the dominant sources 
parameterized by Eq, which would deform the system trajectory away from the depicted solid 
curve, with the deformation becoming small near the fixed points. Conversely, if we wish 
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to follow a randomly chosen trajectory in the phase plane, which is not one of the depicted 
solid lines, we must provide the correct sources in order to ensure that the trajectory is a 
solution of the full gravitational system. But in such cases, the flow lines will be different 
than those depicted, particularly if the new sources turn out to be dominant. 

As a consequence of this discussion, we see immediately that if we want to follow a 
trajectory that would take us to the coordinate origin (/, I) = (0, 0), we must dial the value 
of Eq to zero as time goes on. This of course means that the fluid dominating the universe 
will be neither Hagedorn, nor the momentum mode, nor the winding mode gas. Further, as 
long as the energy density remains positive, the system will have to follow a trajectory which 
is composed piecewise of the infinitesimal portions of the depicted hypebolas, because of the 
last of Eqs. (134T) . This then shows that for such evolution the coordinate origin will be an 
accumulation point of the fixed points on the one-parameter family of hyperbolas enclosed 
between the dashed straight lines in the Figures (jlJ)-(jSJ), and so, it itself will be a fixed point. 
Therefore, it will take an infinite time to reach it. As a result, the universe represented by 
such solutions will be future geodesically complete, and its passage through the origin will 
be impossible. This shows that the idea of the loitering universe dynamics of [14] cannot be 
realized with positive sources. The transition from one side of the phase diagram to another 
will never complete, and the universe on one side, approaching the fixed point, will be left 
hanging there forever. 

One may also ask how robust are these conclusions, since quantum corrections may 
violate positivity of energy, and perhaps relax some of the restrictions on the dynamics. In 
fact one does expect additional contributions to the effective Lagrangian, that could arise 
from other stringy matter, higher order a' or g s corrections, spatial curvature, and various 
possible contributions to the effective dilaton potential, and so this question is relevant. If 
the effective energy density in the course of evolution were negative for a while, then a change 
from (— ) branch solutions to (+) branch ones may occur, redirecting the flow of solutions 
from weak to strong curvature regime [22]. However the transitions from (+) to (— ) branch, 
that could avert the future singularity once the cosmology starts evolving towards it, cannot 
happen if the sources which dominate the universe do not violate NEC during the transition 
from different string phases. To see this, we can generalize Eqs (IT2"|) - (IT4"|) to nominally include 
additional contributions to the dynamics in the effective Lagrangian C m = C m ((f), g^, . . .) 
(denoted by . . .), arising from any other stringy matter, higher order a' and g s corrections, 
spatial curvature, and various possible corrections to the dilaton potential. Restricting the 
analysis to the homogeneous and isotropic perfect fluid cosmologies as before, with the 
string ensemble sources T^ u = diag (p s ,p s , . . . ,p s ), now including the contributions from the 
corrections to C m as well as the string ensemble terms, we find [24] 

yj s = ±y/NH* + e*-p a , (36) 

H s = ±H 8 y/NH* + e<l>°p a + V* (p. + A^C m ) , (37) 

p s = —NH S (p s + p s ) - (fisA^Cm , (38) 

(p s = NH 2 S + l e +* {p s - A £ m ) , (39) 

where A^£ m is the variational derivative of the effective Lagrangian with respect to the 
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dilaton including any couplings of the dilaton to the general sources which are absent in 
the string ensemble description ([3]). Here clearly the upper sign refers to the (+) branch 
configurations, and the lower sign to the (— ) branch. The Eq. (15§)) follows by taking a 
derivative of (136)) and straightforwardly manipulating the terms. 

As noted before, a branch change may not occur unless the discriminant, or the 'egg' 
function [23] 

e = NH? + e*'p a , (40) 

vanishes during some stage of the evolution, which requires p s < 0. Now, suppose that p s 
does dip below zero for some range of variables and parameters. From Eqs. (136)) - (139)) it is 
clear that the phase space trajectories can only explore the region of p s < down to the 
contour e = 0. The region e < is excluded by demanding the reality of the solutions. 
Now, the region where p s < should be compact, in order to match it asymptotically to the 
string gas cosmologies discussed above, for which p s > 0. Thus, the contour e = should 
be bounded, with its interior excluded, just like the 'egg' region discussed in the pre-Big 
Bang attempts to induce branch changing [23,24]. The mechanism of branch- changing can 
then be viewed as a collision between the system trajectory in the phase space and the 'egg' 
contour e = 0, that needs to stand in its way. A branch changing would then require the 
following ingredients: 

• Collision with the 'egg'; branch changing without collisions is impossible, since the sign 
of (p s will not change unless e = somewhere along the trajectory. 

• At the moment of collision, ip s oc ^/e = 0. Since the trajectory cannot enter inside the 
'egg' and will not stop there, as (p s ^ 0, it will ricochet away, grazing the 'egg'. By 
continuity, the sign of ip a , and consequently the branch, will change. 

• Subsequently the trajectory should fly away from the 'egg', such that further collisions 
that can change the branch again will not occur. Alternatively, for the branch change 
to occur there could be only an odd number of collisions. 

So to have a branch changing, we must generate an 'egg' region in the way of a trajectory 
which asymptotically approaches those of Figures OS])-©. Then we must tune the initial 
conditions such that a trajectory will glance off the 'egg', and make sure that it can flow 
away after an odd number of hits. To check when these requirements may be realized, write 
the first derivative of <p s = ±y/e using (1391) . ±4 (>/e) = e + |e^ s p s — |e^ s A^£ m , and using 
(157)) eliminate A^,£ m to find 

± j t (Vi) = -H. + H s ^ s + V* (p s + p s ) . (41) 

Suppose that a branch change does occur, so that the trajectories arriving to the surface of 
the 'egg' at say the instant t^, when e = 0, and departing away from the region of the 'egg' 
at say t e when p s = and so e = NH^, and beyond which again p s > 0, are on different 
branches. Integrating Eq. (|41)) between these instants yields [22] 

\ f e dte*°(p s +p s ) . (42) 



(It 
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The first term is a total derivative and so is given by the difference of the quantity ±y/e + H s 
at the limits of integration. The second quantity is dt H s <p s = f f e d(j) s H s , and so it is 
given as the line integral of H s (4> s ) over the phase space trajectory. Because the flow of 
the phase space trajectories around the 'egg' is clockwise [22], which can be glimpsed from 
the flow lines on the Figures fllJl-flSD, the relation <f) s = (fi s + NH S , and the fact that around 
the 'egg' H s cannot change the sign since it obeys NH] = e + ef\p s \ > ef\p s \ > 0, this 
integral is equal to the area A between the projection of the phase space trajectory onto 
the (<j) s ,H s ) hyperplane and hence it is positive definite. Then integrating and using that 
e(t e ) = NH 2 s (t e ) and e(t h ) = 0, we get 

(±VN-l)H.(t e )+H B (t h )+A=-~ r dte+'ips+ps) . (43) 

Note that this formula immediately applies for any odd number of hits between t% and t e , 
thanks to the additivity of definite integrals. 

Now, if the incoming trajectory is a (— ) branch expanding universe with H s > 0, keeping 
the lower sign in the first term of Eq. (1431) we can rewrite this equation as 

(VN + l) H,(t e ) = H s (t h ) +A+^ r dte*°(p s + p.) , (44) 

which in principle allows a branch change (— ) — > (+) to occur even when the NEC is always 
satisfied, p s + p s > 0. For such a transition, what is required is that prior to the transition, 
the expansion of the universe on the (— ) branch, as characterized by H s (th), is slowed down. 
This is accomplished by the effective string frame energy density p s dipping below zero, but 
then NEC need not be violated. Such transitions can be realized by even such mundane terms 
as the positive spatial curvature of the universe, as is familiar in the examples considered 
in [22]. However, such transitions are not desired. On the contrary, if a transition from (— ) 
branch to a (+) branch occurs, it then implies that the universe will end up in the future 
singularity, unless another transition back to (— ) branch happens. 

Those transitions, which are also required in order to realize the cosmological scenarios 
discussed in [5, 14], are much harder to accomplish. Indeed, going back to Eq. (J4"3"j) . to 
describe the (+) — ► (— ) transition we now must keep the upper sign in the first term, as 
that describes an initially (+) branch solution. Then we can rewrite the resulting equation 
for an expanding universe as 

(VN - l) H s (t e ) + H s (t h ) +A=~ I" dte+-(p. + p s ) . (45) 

Since the universe should be expanding before and after the transition, the left hand side of 
this equation is positive definite, and unless NEC is violated, we see that no solutions can 
satisfy it. Thus, if any solutions that are consistent with Eq. (j4"5l) were to exist, NEC must 
be violated in an overwhelming way - during the transition, the universe must be dominated 
by the sources which do not respect NEC. 

It is very difficult to imagine such sources emerging in string theory. It is true that 
Casimir energy violates NEC, however in cosmological applications such contributions are 
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by and large sub leadings to the dominant sources that control the evolution [36]. One 
exception, where the NEC violating effects are studied and found to significantly affect the 
background is the eternal inflation [38]. In eternal inflation, quantum fluctuations drive the 
scalar field up the potential, which violates NEC, and induces an inflationary landscape which 
populates the vacua of the theory. This happens without violating the positivity of energy 
(while keeping the quantum matter of the universe in equilibrium with gravity, as opposed 
to just in equilibrium) in contrast to what branch changing processes would need to do. On 
the other hand, the recent explorations of the string landscape have shown that it is possible 
to find inflation in string theory [39], and so such dynamics may occur even if the universe 
starts out dominated by some string gas. If this happens, however, inflation will generically 
completely wipe out the traces of the early initial conditions in the universe, reshaping the 
universe in accord with inflationary dynamics. This would then solve cosmological problems 
without the need to resort to the gas phase. 

A possible method to induce NEC violation is to introduce ghosts in the theory, which 
of course typically lead to a host of problems with both the microscopic and the large scale 
behavior of physical systems, involving uncontrollable instabilities. An interesting suggestion 
that a ghost may condense [26], and its instabilities tamed in cosmological applications 
has been made in [27,28], however string theory has been proven to be ghost-free [29,30]. 
Furthermore, the ghost condensate model [26] is known to have a very low cutoff, which 
is necessary to keep the bad behavior under control. To go beyond the cutoff and explore 
energy scales such as those encountered in string cosmology, one should need to find a UV 
completion of the theory and use it at the scales above the cutoff. The investigation of the 
possible completions however has turned up arguments against such completions based on 
locality and causality [31]. Thus unless it is demonstrated that the ghost condensate can 
be consistently embedded in string theory, the interesting condensed ghost bounces of [27] 
do not really provide the mechanism to generate branch changes in string gas cosmologies. 
Thus for now, it appears that the phenomenological applications of (+) branch solutions, 
that evolve towards future singularities, invariably leads one into the string swampland [40] . 

4 Summary 

To summarize, in this paper we have considered the obstructions to phase transitions between 
different thermodynamic phases of the string gas in the early universe. These arise because 
the string dilaton and gravity react to the background sources by generating the flow of the 
dilaton, or the string coupling. If the energy density of the sources is non-negative, the sign of 
the time derivative of the dimensionally reduced dilaton (p s = <p s —N\ s is conserved, behaving 
like a discrete charge, and separating the solutions into two superselection sectors, or (+) 
and (— ) branches. The (+) branch solutions are problematic since they evolve to a future 
spacelike singularity, which looks either like a Big Crunch or a pole inflation in the string 
frame. The singularity separates different branches of string gas solutions like a geometric 
precipice in the moduli space, and prevents them from joining together. The flow cannot 

2 The corrections to the Hagedorn gas equation of state which have been analyzed in [37] indeed do not 
violate NEC. 
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be altered without negative energy and NEC violations. This excludes the early universe 
cosmologies based on the transitions of the string Hagedorn gas evolving away from the 
fixed point to the FRW radiation cosmology [5], and the loitering gas models of [14], unless 
new sources are introduced which violate both positivity of energy and NEC. We provide 
quantitative criteria which such sources must satisfy to yield a branch change, but recall 
that the known mechanisms which could in principle satisfy these criteria [27, 28] do not at 
present come out of string theory in a natural way. It does remain possible to build composite 
cosmologies using the solutions that reside on the same branch, and we have seen an example 
of a model patched together from winding mode, Hagedorn and momentum mode cosmologies 
on the (— ) branch. Those solutions behave differently from the ones discussed in [5, 19] 
because they do not have a future singularity. The past singularity could be 'resolved' 
by using T-duality, instead of bounces induced by negative energy and NEC violations, 
and is closer in spirit to the original ideas of how T-duality might help in cosmological 
applications. These solutions do not yet have interesting cosmological applications, but it 
would be interesting to consider what happens with cosmological perturbations in universes 
that are T-dual to each other. 
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